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Abstract 

We prove that the stochastic differential equation 

ft — s 

(f) Y 3 , t (x) = Y 3 , 3 (x)+ / f(Y s ,s+u(x)) dX s +u, 
Jo 

driven by a Levy process whose paths have finite p-variation almost surely for some 
p £ [1,2) defines a flow of locally C 1 -diffeomorphisms provided the vector field / is 
a-Lipschitz for some a > p. Using a path-wise approach we relax the smoothness 
condition normally required for a class of discontinuous semi-martingales. 



Introduction 

In this paper we give sufficient conditions for the solution of the following in- 
tegral equation to define a flow of local C 1 -diffeomorphisms on M. d when the 
integrator is allowed to have unbounded variation: 

(2) Y Stt {x) = Y s<s {x) + / f{Y StS+u (x)) dX s+u , Y StS {x) =xe R d . 



The above assumption means that one must take care defining the integral. 
In this paper we consider integrators with finite p- variation for some p£ [1,2). 
We use the Young JlTj integral 

(3) / f(u)dg(u) (Young), 



o 



which is defined by taking a limit of Ricmann sums whenever the functions / 
and g have no common discontinuities and have finite p and q variation for some 
l/p+ l/q > 1. Lyons solved the following integral equation: 



(4) Vt = yo+ / f(y u ) dx u y = aeM. d 



o 

where Xt is a continuous function, having finite p-variation for some p < 2 and 
/ is a vector field which is a-Lipschitz for some a > p — 1. Furthermore, he 
proved that if a > p then the solution is unique. 
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In this paper we show that the solutions of equation (Q) form a flow of local 
C 1 -diffeomorphisms. Then we extend the solution of equation (Q) to allow the 
integrator to be discontinuous. Finally we apply the results to SDEs driven by 
a class of Levy processes whose paths, almost surely, have finite p- variation for 
some p < 2. It is shown (Corollary |4.l| ) that one can relax the smoothness con- 
dition used when using the ltd calculus within the semi-martingale framework. 

The paper is organised in the following way: The first section introduces 
the notation used in the remainder of the paper. The second section contains 
the proof of the diffeomorphic flow property. In the third section we extend the 
result of Lyons to incorporate integrators with discontinuities. Finally we apply 
the results to the class of Levy processes whose paths a.s. have finite p- variation 
for some p < 2. 

1 Preliminaries 

In this section we recall some definitions which will be required throughout the 
paper. 

Definition 1.1 The p-variation of a function X(s) over the interval [0, t] is 
defined as follows: 

(5) ll*IU.*] = { sup ^2\X(t k )-X(t k -x) 

where 7r[0, t] is the collection of all finite partitions of the interval [0,t]. 



Remarks 1.1 Note that \\.\\ p is a semi-norm on the class of functions whose 
p-variation is finite. Let W p denote this class of functions. Trivially one has W p 
C W q when p < q. Let W p (a) denote the subspace of paths in W p with initial 
point a. Then W p (a) equipped with \\.\\ p is a Banach space. 

Let X t be a path of finite p-variation on [0,T]. Then necessarily there are 
only countably many discontinuities and there exist left and right limits of X t 
for all t S [0, T] . 

Note that we use a stronger version of p-variation than is used in the stochas- 
tic calculus literature. Fristedt compares the different definitions of variation. 

Example 1.1 Probabilists are used to the idea that Brownian motion a.s. has 
finite quadratic variation. The process a.s. fails to have finite 2-variation in the 
above (strong) sense; it does however have finite p-variation a.s. for all p > 2. 

Definition 1.2 A function f is in Lip (a) for some a > 1 if 

df 

II /II oo < oo and ^— G Lip(a - 1) j = 1, . . . ,d. 
Its norm is given by 

d df 

||/||Lip(a) = II/IU + 5Z WqT. Lp(a-1) f0Ta > h 

j=l 3 
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This is Stein's |T^] definition of a-Lipschitz continuity for a > 1. It extends 
the classical definition: / is in Lip(a) for some a £ (0, 1] if 

\f(x)-f(y)\<K\x-y\ a 

with norm 

l/l. t ,M. 

x ^y \x-y\ a 

2 Diffeomorphic flows 

Consider the following differential equation: 

d 

(6) dy t =Y J f\yt)dx\ y = aeR d 

i=l 

where (xt)t>o is a continuous R"-valued function with finite p- variation for some 
p < 2 and / : R d — > hom(M™, Lyons proved the following theorem: 

Theorem 2.1 (@) J// is a-Lipschitz for some a > p — 1 f/ien solutions to 
(||) exisi. J/ / is a-Lipschitz for some a > p then there is a unique solution. 



Remark 2.1 The existence result proves that the solution path has finite pf- 
variation for some pi > p. The uniqueness result proves that the solution has 
finite p-variation. Both results are proved by applying fixed point theorems to a 
sequence of Picard iterations. 

We modify the proof of the above theorem to show that a flow of local 
C^diffeomorphisms exist. First, we recall the definition of a flow of diffeomor- 
phisms. 

Definition 2.1 A collection of maps {<j> s ,t '■ < s < i] form a flow of homeo- 
morphisms on M. d if: 

(i) 4>s,u = 4>s,t o 4>t,u for all s <t < u. 
(ii) 4> StS is the identity map on R d for all s £ R + . 
(Hi) <p St t is an onto homeomorphism ofW 1 for all s,t £ K + . 

Further if (f> s .t satisfies (iv) then it is called a flow of C k -diffeomorphisms. 

(iv) (j) Sy t{x) is k times differ entiable with respect to x. The derivatives are con- 
tinuous with respect to (s,t,x). 



Remark 2.2 Note that a stochastic flow of homeomorphisms ^ P-H4] * s a 
continuous random field which satisfies the above properties off a set of measure 
zero. 
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Theorem 2.2 Let (X t )t>o be a M. n -valued continuous path having finite p- 
variation for some p < 2. Let f : M. d — > hom(R™, K d ) fee a-Lipschitz for some 
a > p. Then the unique solution to the following system of differential equations 
is a flow of local C 1 -diffeomorphisms on M. d : 

ft— 8 

(7) Y Sit (x)=Y SiS (x)+ / f(Y s . s+u (x)) dX s+u , y s , s (i)=i£l d . 

Jo 

The proof is an extension of the method used in jlo| . The key step is to consider 
the following pair of differential equations: 

dY s . u (x) = f(Y s , s+u (x)) dX s+u Y s>s (x) =xeR d 

(8) dK s , u (x) = Vf{Y s , u (x)) K s . u {x) dX s+u K s>s {x) = Ie M dxd 

where M dxd denotes the space of d x d matrices over M. 
Proof. We rewrite (||) as follows: 

(9) dZ s>u (x) = h(Z s , u {x)) dX s+u Z StS (x) = {x, L) £ R d x M dxd , 
where h is the vector field defined by 

(10) h((y,k))±(f(y),Vf(y)k). 

Given that / is a-Lipschitz for some a > p we know that h is (a — 1)-Lipschitz. 
By applying Lyons' existence result jl0| we deduce that there exist solutions to 
([)]). However, the solution of the first equation in (||) is unique. This implies 
that the solution of (||) is unique. 

It remains to show that the solution is C . Let a, b £ M. d . Let q £ (p, 2) and 
r = (a — l)/p — S where 5 £ (0, (a — p)/p). Then 

\\K s , t (a) - K Stt {b)\\ p ^ T] < 
2 (p-i)/J|| /"(V/(F S)S+M (a)) - Vf(Y 8 , s+u (b)))K s , s+u (a) dX s+u \\ p ^ T] 



+ || / Vf{Y 3 , s+u (b))(K s , s+u (a) - K SjS+u (b)) dX s 

+u\\p,[0,T] (■ 



Let 

{K StS+u (a) - K S}S+u (b)) dX s+v 



h= K StS+u (a) dX s+u , g t = 
Jo Jo 

A bound of Young implies that 

IHI P ,[0,T] < (l + ((l/p+l/q))\\K s , s+ Xa)\\ q ^ T] \\X\\ p , [0iT] , 

\\9\\ P ,lo,T] < (l + C(l/p + l/9))||^.,.+.(a)-^a,.+.(6)||,,[o,2l IWIp,[0,T]- 

So 

\\K s<t (a) - K s<t (b)\\ pAOjT] < 

2(p- 1 )/f||X|| Pi[0iT] (i + C(i/p + i/g)) 

x |||V/(F S)S+M (a)) - V/(r,, a+tt (6))|| r)[0l ril|iifa,.+.(a)ll« > [o J ri 

+ ||V/(y s , s +„(6))||,,[o,T]||^ s , s +.(a)-K s , s+ .(&)|| 9 ,[o,T]] 
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Noting that 
\\Vf(Y StS+u (a)) - V/(F 5)S+ „(6))|| g , [0 ,T] 
< 2|lV/Cn,-+»(a)) - V/(n,.+™(6))H^ [ S, ri ll/ll&r )/fl! II* r («) - * r (6)llS|- ,, * )/9 
where (?* G (p,q), we see that for T sufficiently small 

\\K s<t (a) - K a!t (b)\\ p , [0! T] 

varies continuously in space. 

We have not yet considered the inverse maps and their regularity. The above 
proof can be modified to recover the result by considering the solution of the 
differential equation driven by the time-reversed path (X t - u )o<u<t- D 

We note that the pathwise approach improves the following result of Kunita 
HI whenever the continuous integrator almost surely has finite p-variation for 
some p < 2: 

Theorem 2.3 ([^] Theorem 2.3) Let f be a C k map with a bounded deriva- 
tive, let (x t )t>o be a square integrable martingale. Then the solution of (||) 
defines a C k ~ 1 flow of diffeomorphisms for all t with probability one. 

The improvement is due to the increased smoothness of the integrator. It is 
the pathwise approach that picks up the p-variation (p < 2) of the integrator 
thereby reducing the smoothness requirement on /. 

The following example contains an interesting application of the above result. 
The example, fractional Brownian motion, is a process which is not a semi- 
martingale H Corollary 2.2]. The process had been considered suitable for 
financial modeling due its heavy tail distribution. (However Rogers |L4| has 
shown that arbitrage opportunities appear when using the fBm.) 

Example 2.1 Let (X t )t>o be a fractional Brownian motion with Hurst index 
H e (i,l). That is, Xt is a mean zero, Gaussian process whose covariance 
function is given by 

(11) T H ( S , t) = \{t™ + s 2H \t sf H ) Var{ Xl ) 
where 

(12) VariX,) = — + J ((1 + s) H ^ s^f ds. 

Lemma 2. 1 Fractional Brownian motion with Hurst index H has finite p- 
variation almost surely for p > 1/H. 

Proof. We estimate \X t — X s \ p for s < t £ R + as follows: Partition the 
interval [s,t] by selecting the largest dyadic interval [{k — l)2 _n T, k2~ n ] which 
is contained within [s,t]. Then we add dyadic intervals to either side of the 
initial interval, which are chosen maximally with respect to inclusion in the 
interval [s,i\. Continuing in this fashion we label the partition according to 
the lengths of the dyadics. We note that there are at most two dyadics of the 
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same length in the partition we can label these [Zi,fc , '"l.fc] and [h,k, r 2,k] where 
n,k < h,k ■ Then 

oo 

(13) [u,v] - (J U ikk,n, k }. 

fe=l i=l,2 

We use the Holder inequality to show that 



(14) \X t -X.\*< E E \ X h,*- X nJ 

\ fc=l i=l,2 

/ OO s p— 1 OO 

^n=l ' n=l ^i=l : 2 

oo 

(16) <C(p ) T)E n7 El I '«,*- I n,J P - 

n=l i=l,2 

One can uniformly bound \X t — X s \ p for any pair of times s < t g [0, T] by 
extending the above estimate over all dyadic intervals at level n, that is 

oo 2" 

(17) \X t -X s \P<C(p, 7 ) ^nf^X^-Xr^f. 

n—l i—1 

The p- variation of X t can be estimated by this bound. 

oo 2™ 

(is) sup £ |x tfe+1 - x tfc r < c(p, 7) E n7 E - i p - 

7r n—l i—1 

Noting that 

\2H 



(19) E^(X t - X s ) z j =\t- s\ ZM VariXx). 

we see that the right hand side of ( |18|) is finite a.s. whenever Hp > 1. □ 
From Theorem 2.2 and Lemma |2.1| we deduce the following result: 

Corollary 2.1 Let (X t )t>a be a fractional Brownian motion with Hurst index 
H. Let f be an a-Lipschitz function for some a > l/H. Then the solution of 
(0) driven by (X t ) defines a stochastic flow of local C 1 -diffeomorphisms. 



3 Discontinuous paths 

In this section we extend the results of Section || to discontinuous integrators. 
The theory is applied to sample paths of certain classes of Levy processes, 
namely those that have finite p-variation a.s. for some p < 2, see Section ^. 

First we determine the solution's behaviour when the integrator jumps. 
There are two possibilities which we consider: the first is an extension of the 
Lebesgue-Stieltjes integral, the second is based on a geometric approach. 

Suppose that the discontinuous integrator has bounded variation. Then the 
solution Y would have a jump of the form 



Y t -Y t -=f{Y t -) {Xt-Xt-) 
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at a jump time of X . If X has finite p- variation for some 1 < p < 2 we could 
insert the (countable) set of jumps of the above form at the discontinuities of 
X. The disadvantage of this approach is that one cannot expect the solution to 
generate a flow of diffcomorphisms for a general / 

Another jump behaviour which we use is the following: When a jump of the 
integrator occurs we insert fictitious time, allowing the jump to be traversed, 
thereby creating a continuous path. Then we solve the differential equation 
driven by the continuous path. Finally, we remove the fictitious time compo- 
nents of the integrator and the solution path. We call this a geometric solution. 
This jump behaviour has been considered before by (ll[] and H. The reason 
that we work with this second type of solution is that it experiences an 'in- 
stantaneous flow' along an integral curve at jump times. This removes the flow 



problems associated to the discontinuous integrator 12j , || . 

Now we define a parametrisation for a cadlag path X of finite p-variation. 
The parametrisation adds fictitious time, allowing the traversal of the discon- 
tinuities of X. We prove that the resulting continuous path has the same p- 
variation as the original path. We solve (||) driven by the continuous path 
from which we extract a geometric solution by removing the fictitious time (or 
undoing the parametrisation). 

Definition 3.1 Let X be a cadlag path of finite p-variation. Let 8 > , for 
each n > 1 , let t n be the time of the n 'th largest jump of X . We define a map 
t s : [0,T] — > [0,T + S^iLi \ J{ti)\ p ] (where J{u) denotes the jump of the path 
X at time u) in the following way: 

oo 

(20) r\t)=t + SY / \J(tn)\ p X{t n < t} (t). 

n=i 

The map t 5 : [0,T] -> [0,t 6 (T)} extends the original time frame into one on 
which we define the continuous process X s (s) below. 

Definition 3.2 



x s (s) 
(21) = 



x(t) ifs = r s (t), 

<K) + (s T S (t-))j(t n )5-i\j(t n )\-P ifs G [r s (t7 n )r s (t n )). 



Remarks 3.1 We see that (s, X*), s e \0,t 5 (T)] is a parametrisation of the 
driving path X . The terms \ J(t n )\ p in ( J20| ) ensure that T 5 (t) remains finite. 



The following two results are proved in |lq ] : 

Proposition 3.1 J7^ , Proposition 1.1] Let X be a cadlag path of finite p- 
variation. Define a parametrisation of X as above, then 



( 22 ) 11^ \\p,[o,r s (T)] = ll^llp,[o,T] VJ > 0. 

Theorem 3.1 J7^ , Theorem 1.2] Let X be a cadlag path with finite p-variation 
for some p < 2 . Let f be a Lip(7) vector field on W 1 for some 7 > p . Then 
there exists a unique geometric solution Y with finite p-variation which solves 
the differential equation 

(23) dY t = f(Y t ) dX t Y Q =ae R". 
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We use Theorem 2.2 to deduce 

Corollary 3.1 Let X be a cadldg path with finite p-variation for some p < 2 . 
Let f be a Lip(7) vector field on W l for some 7 > p . Then the unique geometric 
solution to the following system of differential equations is a flow of local C 1 - 
diffeomorphisms on M. d : 

(24) Y s<t (x)=Y S!S (x)+ f(Y s , s+u (x))dX s+u , Y s , s (x) = x e R d . 



4 p-variation of Levy processes 

In this section we recall some results concerning the p-variation of Levy pro- 
cesses. These determine when Corollary 3.1 can be applied to sample paths of 
Levy processes. 

Levy processes are a combination of drift, diffusion and jumps. The drift 
component has bounded variation whilst the Gaussian part has finite p-variation 
a.s. only for p > 2. Therefore we restrict our attention to the compensated, 
pure jump process given by: 



(25) lim < / [x N t {dx) - tx v(dx) 

£ ->° I Je<\x\<l \ 

To characterise the p-variation of Levy processes it is natural to look at the 
jump behaviour. Blumcnthal and Getoor measured this in the following way: 

Definition 4.1 ||/ The index of a Levy process X is defined to be 



(26) /3 = hrNa>0: / \y\ a v{dy) < 00 

I J\y\<l 

It is a measure of the frequency of the process' small jumps. If the integral 

(27) / \x\ v{dx) 



is finite then we see that the compensated jump process has bounded variation. 
In |l[, Blumenthal and Getoor proved the following result for the symmetric 
a-stable processes a <E (0, 2] : 

Theorem 4.1 Theorem 4-1] Let {X t ; t > 0} be the symmetric stable process 
in R" of index a G (0, 2] . Then 



(28) PM|X||^ [0> t] = ooJ =1 or 

depending on whether [3 < a or [3 > a . 

Their proof centred on the special form that the characteristic functions of a- 
stable Levy processes (a £ (0,2]) take: 

iuXt 1 —\tu a 



(29) E 

It was conjectured that any pure jump Levy process would have finite p- 
variation provided p > (3. Bretagnolle |3| resolved the issue when the index 
(3 < 2 with the following result: 
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Theorem 4.2 ffl LetX(t) be a process with independent increments with Levy 
measure v . Then X has finite p-variation for some p < 2 on [0, 1] a.s. if and 
only if X has no Gaussian part and 

(30) f \x\ p v{dx) < oo. 

J\x\<l 



Remark 4.1 In \13jl the characterisation of the sample path p-variation of all 
Levy processes is completed. 



Using Theorem 4.2 we apply Corollary 3.1 to a class of Levy processes. 

Corollary 4.1 Let (X t ) be a Levy process without a Gaussian part and with 
index (3 less than two. Let f be an a-Lipschitz vector field for some a > ft. 
Then there exists a unique stochastic flow of local C 1 -diffeomorphisms (Y s , t ) 
which solve the equation: 



(31) Y S;t (x) = Y s>s {x) + 



f(Y s<s+u (x)) dX s+u , Y s>s (x) = x. 
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